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Supermassive black holes with masses of millions to billions of solar masses are commonly found in the centers of 
galaxies. Astronomers seek to image jet formation using radio interferometry but still suffer from insufficient angular 
resolution. An alternative method to resolve small structures is to measure the time variability of their emission. Here 
we report on gamma-ray observations of the radio galaxy IC 310 obtained with the MAGIC telescopes, revealing 
variability with doubling time scales faster than 4.8 min. Causality constrains the size of the emission region to be 
smaller than 20% of the gravitational radius of its central black hole. We suggest that the emission is associated with 
pulsar-like particle acceleration by the electric field across a magnetospheric gap at the base of the radio jet. 
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More than three decades ago it was proposed that the radio emission of 
extragalactic jets results from a relativistically moving plasma consisting 
of magnetic fields and accelerated particles following a power-law energy 
distribution (7). One of the major assets of the model is that it can explain 
the non-thermal emission of extragalactic jets across the entire electro- 



magnetic spectrum, from radio waves up to gamma rays. The emission 
can be understood as synchrotron radiation and inverse Compton scat- 
tering (2, 3) due to particles accelerated at shock waves in the jets. The 
gamma rays can reach very high energies measured in Giga-electronvolts 
( 1 GeV =10^ eV corresponding roughly to the rest mass energy 
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equivalent of the proton) and Tera-electronvolts (1 TeV= 10^^ eV). 
According to the Blandford-Znajek mechanism, the jets are powered by 
extracting rotational energy from the black holes which have acquired 
angular momentum through the accretion of surrounding gas and black 
hole mergers (4), although so far astrophysical evidence for the role of 
black hole spin is still lacking (5). For a maximally rotating supermassive 

black hole of mass , where MO denotes one solar 

mass, the size of the jet formation region should be of the order of its 

gravitational radius = G^M / 1 .5 X 10^^ ffl^ m and twice 

this value for a non-rotating Schwarzschild black hole. Astronomical 
telescopes do not yet provide the angular resolution needed to image 
structures on this scale. The highest resolution images of jets obtained 
with very long baseline radio interferometry show radio-emitting knots 
traveling down the jets (6). Approaching the black hole, the spectra cut off 
at increasingly higher frequencies due to synchrotron-self-absorption. 
Observations at very high frequencies where the core becomes transparent 
are needed to zoom into the region where the jets are emerging from. The 
record holder is a very long baseline radio interferometry observation of 
the jet of the nearby radio galaxy M87 at a frequency of 230 GHz, re- 
solving a source with the size of 1 1.0 i 0.8 gravitational radii (7). 

The event horizon light-crossing time 

While direct imaging of the jet formation region has to await better 
angular resolution, indirect information about its size can be inferred from 
the temporal variability of the emission coming from that region. The 
observed gamma-ray variability time scales indeed reach down to the 
event horizon light crossing time 

^^BH = / <^ = I =8.3 minutes {8) vindicating the 

scenario that the jets originate from the magnetospheres of accreting black 
holes. An example is the radio galaxy M87 in the Virgo cluster of galaxies 
(P, 10). This galaxy harbors a central supermassive black hole with the 

enormous mass of-^ 6.4x10 Mo (7). M87 exhibits gamma-ay 
variability at a time scale of days (72) which is consistent with the 

light-crossing time of the event horizon A^gj^ (M87) = 0.4 days. 
Smaller than a black hole? 

Observations of Mrk 501 {13) and PKS 2155 - 304 {14) at very high 
energies have provided evidence for extreme variability events with flux 
doubling time scales as short as ~ 2 min. The ultrafast variability cor- 
responds to 1st < A^gjj and therefore casts a shadow of doubt on the 

current shock-in-jet paradigm. It has been suggested that relativistic bulk 
motion of the jets could explain the observed time scales (75). The ar- 
gument relies on the observation that these flaring sources belong to the 
class of sources that astronomers call blazars. In blazars, the jets are 
pointing at a small angle toward the observer. Since the jet plasma moves 
with a speed close to the speed of light (3 = v/c ^ 1 leading to a bulk 

Lorentz factor Fj > 1 , several effects arise due to the relativistic 

boosting of the emission (7(5). One of them affects the time scale of flux 
variations of the emission from a shock. The moving shock plasma almost 
catches up with its own radiation, and this leads to a shortening of the 

observed variability time scale /St compared with the variability time 

scale A^' in a frame comoving with the shock given by 

1st = {\ + z)5 ^ ISf where Z denotes the cosmological redshift of 



the source. For a given value of Fj , the Doppler factor 5 depends 

strongly on the orientation angle of the jet 0 (Fig. 1). Note that ^ = 0 
corresponds to perfect alignment. For Mrk 501 and PKS 2155 — 304, 
almost perfectly ahgned jets with Fj > 50 would be needed to ac- 
commodate for 1st < A^gjj and to avoid self-absorption of the gamma 

rays due to pair production (75). In blazars, interferometric observations 
of the superluminal motion of radio knots suggest lower values 

F j ~ 1 0 and orientation angles of a few degrees (7 7). Larger values of 

Fj would lead to a problem with population statistics: the number of 

unbeambed counterparts of blazars viewed at larger angles would then 
exceed the number of radio galaxies, commonly believed to represent the 
misaligned blazars {18). Assuming lower black hole masses would bring 

down A^gjj . However, lower masses conflict with the firmly established 

dynamical measurements of black hole masses, and therefore do not seem 
to be a likely solution of the dilemma. Other possible solutions of this 
Doppler factor crisis {19) invoke models of structured jets {20, 21) or 
Poynting flux dominated jets in which only a few but very fast seed par- 
ticles at the jet base reach high Lorentz factors {22-24), before the Poyn- 
ting flux is converted into the kinetic energy of the bulk flow by mass 
entrainment. 

All of these attempts to explain the sub-horizon scale variability with 
relativistic projection effects alone ignore a fundamental problem. If the 
perturbations giving rise to the blazar variability are injected at the jet 
base, the time scale of the flux variations in the frame comoving with the 

jet is affected by time dilation with Lorentz factor Fj . In blazars where 
^5^ ~ F j , the Lorentz factor cancels out, and the observed variability time 
scale is ultimately bounded below by A^gj^ . 

IC 310: A gamma-ray lightning inferno 

IC 3 10 is a peculiar radio galaxy located in the outskirts of the Perseus 
cluster at a distance of 260 million light years from Earth. On November 
12/1 3th in 2012, MAGIC, a system of two Imaging Atmospheric Che- 
renkov Telescopes located on the Canary island of La Palma (25), de- 
tected an extraordinary outburst of gamma rays from this object (see Fig. 
2, SI). The details of the analysis can be found in section S2. Prior to these 
observations, variable gamma-ray emission from IC 3 10 had already been 
detected by satellite and ground-based gamma-ray instruments at GeV 
and TeV energies, e.g., Fermi-LAT and MAGIC {26-28). In the night of 
the flare in November 2012, the mean flux above 300 GeV was 

(6.08 ± 0.29) X 10 " cm s , i.e., four times higher than the 

highest flux during previous observations in 2009/2010. The measured 
spectrum (Fig. 3) can be described by a simple power law with a differ- 
ential photon spectral index of F = 1 .90 ± 0.04^^^^ ±0.1 5^^^^ in the 

energy range of 70 GeV - 8.3 TeV (see Table S2). Owing to its proximity, 
the spectrum of IC 310 is only marginally affected by photon-photon 
absorption in collisions with the extragalactic background light (EBL). 
IC 310 harbors a supermassive black hole with a mass of 

M = (3^)xlO^MO (see section Sl.l) corresponding to an event 

horizon light-crossing time of At^^ =(23^^^^^ min. The mass has 
been inferred from the correlation of black hole masses with the central 
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velocity dispersion of their surrounding galaxies (29, 30). The reported 
errors are dominated by the intrinsic scatter of the distribution. The same 
value of the mass is obtained from the fundamental plane of black hole 
activity (31). The scatter in the fundamental plane for a single measure- 
ment is larger and corresponds to a factor of 7.5 . 

During 3.7 hours of observations, extreme variability with multiple 
individual flares has been detected (Fig. 4, S3, S4). The flare has shown 
the most rapid flux variations ever observed in extragalactic objects, 
comparable only to those seen in Mrk 501 and PKS 2155 — 304. A con- 
servative estimate of the shortest variability time scale in the frame of IC 

310 yields A^/(l + z) = 4.8 min. It is the largest doubling time 
scale with which the rapidly rising part of the flare can be fitted with a 
probability > 5% (see fig. S4). The light curve also shows pronounced 
large-amplitude flickering characterized by doubling time scales down to 
min. The conservative variability time scale corresponds to 

20% of the light travel time across the event horizon, or still 60% of 
it giving allowance for the scatter in the dynamical black hole mass 
measurement. 

From the absence of a counter radio jet and the requirement that the 
proper jet length does not exceed the maximum of the distribution of jet 
lengths in radio galaxies, the orientation angle was found to be in the 

range ^ ~ 10 -20 (see section SI. 2), and the Doppler factor con- 
sistent with S ^ 4 (32). These values put IC 310 at the borderline 
between radio galaxies and blazars. The jet power estimated from ob- 
servations of the large-scale radio jet is L- = 2x1 0^^ erg s ^ as- 
suming that it contains only electrons, positrons, and magnetic fields in 
equipartition of their energy densities (see section SI. 3). For a radiative 

efficiency of 1 0% , the Doppler-boosted average luminosity of the jet 
emission amounts to erg s which is close to 

the one observed in very-high energy gamma rays. For S ^ 4 , the 
variability time scale in the comoving frame of the jet, where it should be 
larger than FjA^gj^ , is actually close to A^gj^ (cf Figure 1). A very 

high value of the Doppler factor is required to avoid the absorption of the 
gamma rays due to interactions with low-energy synchrotron photons, 
inevitably co-produced with the gamma rays in the shock-in-jet scenario. 
The optical depth to pair creation by the gamma rays can be approximated 
by 

(10 TeV) ~ 300( J / 4)'" {At 1 1 min)"' [l^^^ / lO^'erg s"' ) • 

Adopting a non-thermal infrared luminosity of 1 % of the gamma-ray 
luminosity during the flare, the emission region would be transparent to 
the emission of 1 0 TeV gamma rays only if (5^10. For the range of 
orientation angles inferred from radio observations, the Doppler factor is 
constrained to a value of (5 < 6 (Fig. 1). One can speculate if the inner 
jet, corresponding to the unresolved radio core, bends into a just-right 
orientation angle to produce the needed high value of the Doppler factor 
(see section SI. 2). The probability for such an alignment seems to be 
rather low. Moreover, the observed radio jet does not show any signs of a 
perturbation of its flow direction on the parsec and kiloparsec scales. 
Since perturbations of the flow direction of the inner jet would later 
propagate to larger scales, major bends apparently never occurred in the 
past. 

In summary, trying to interpret the data in the frame of the 



shock-in-jet model meets difficulties. Considering the role of time dila- 
tion renders a solution of this problem impossible for any value of Fj . 

Therefore, we conclude that the observations indicate a sub-horizon scale 
emission region of a different nature. 

Possible origins of sub-horizon scale variability 

Substructures smaller than the event horizon scale emitting highly 
anisotropic radiation (to avoid pair absorption) seem to be responsible for 
the minute-scale flux variations. The possible explanations fall into three 
categories: (i) mini-jet structures within the jets (33), (ii) jet-cloud inter- 
actions where the clouds may originate from stellar winds (34-36), and 
(iii) magnetospheric models (37-41), similar to those known from pulsar 
theory. 

According to the mini-jet model (i), plasmoids resulting from mag- 
netic reconnection traveling down the jet with a relativistic speed are 
responsible for the minute-scale flux variations observed in blazars. The 
model could help to mitigate the constraints on the bulk Lorentz factor by 
introducing a larger effective bulk Lorentz factor for the plasmoids. The 
mechanism also predicts reconnection events from regions outside of the 

beaming cone T ■ ^ that could explain the day-scale flares from the 

radio galaxy M87 invoking external radiation fields as a target for inverse 
Compton scattering (33). However, the off-axis mini-jet luminosity de- 
pends on (r j^) ^ and the jet power required for IC 3 1 0 is two orders of 

magnitude higher compared to the one estimated from radio observations 
(see section S3. 2). Thus, this model is challenged by the observed high 
luminosity in IC 3 10 during the flares. 

Substructures smaller than the jet radius may also be introduced by 
considering interactions between clouds and the jet (ii). The original 
shock- in-jet model (1) considered this to be the main source of mass 
entrainment and predicted variability from the process. Recently, more 
elaborate work on the model has had some success in explaining the 
variability of M87 by pp-collisions due to the bombardment of clouds 
boiled off red giants with protons in the jet (36). However, the model is 
linked to the cloud crossing time of the jet and the proton-proton cooling 
time, both of which by far exceed the event horizon scale. Faster varia- 
bility could be observed in case the cloud gets destroyed but a strong 
beaming effect would then be needed to explain the observed luminosi- 
ties. In another variant, drift acceleration of particles along the trailing 
shock behind the stellar wind of a star interacting with the jet is consid- 
ered. This might lead to an extremely anisotropic emission pattern. As 
mass-loosing stars sweep across the jet, passing magnetic field lines 
pointing to the observer, the postulated accelerated particle beams in their 
trails become visible for a short time. For IC 310, the emission would 

have to be confined to within an angle of CC ^10 ^ rad to explain the 
observed variability time scale, requiring a very stable direction of the 
accelerated particle beams, at a large angle to the jet main thrust. Since 
two-fluid particle beams are prone to numerous plasma instabilities, the 
scenario relies on unphysical assumptions. 

In magnetospheric models (iii), particle acceleration is assumed to be 
due to electric fields parallel to the magnetic fields. This mechanism is 
known to operate in the particle-starved magnetospheres of pulsars, but it 
could also operate in the magnetospheres anchored to the ergospheres of 
accreting black holes (see Fig. 5). Electric fields can exist in vacuum gaps 
when the density of charge carriers is too low to warrant their shortcut. 
The critical charge density for the vacuum gaps is the so-called 
Goldreich- Julian charge density. Electron-positron pairs in excess of the 
Goldreich- Julian charge density can be produced thermally by pho- 
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ton-photon collisions in a hot accretion torus or corona surrounding the 
black hole. It has also been suggested, that particles can be injected by the 
reconnection of twisted magnetic loops in the accretion flow (39). A 
depletion of charges from thermal pair production is expected to happen 
when the accretion rate becomes very low. In this late phase of their 
accretion history, supermassive black holes are expected to have spun up 
to maximal rotation. Black holes can sustain a Poynting flux jet by virtue 
of the Blandford-Znajek mechanism (4). Jet collimation takes place rather 
far away from the black hole at the scale of the light cylinder beyond 

~ lOrg . Gaps could be located at various angles with the jet axis cor- 
responding to the polar and outer gaps in pulsar magnetospheres leading 
to fan beams at rather large angles with the jet axis. The gap emission 
must be highly variable, since gap height and seed particle content depend 
sensitively on plasma turbulence and accretion rate. For an accretion rate 

of /iz 1 0 ^ (in units of the Eddington accretion rate) and maximal 
black hole rotation, the gap height in IC 310 is expected to be 
h-0.2r^ (40) which is in line with the observations. Depending on the 

electron temperature and geometry of the radiatively inefficient accretion 
flow, its thermal cyclotron luminosity can be low enough to warrant the 
absence of pair creation attenuation in the spectrum of gamma rays. In this 
picture, the intermittent variability witnessed in IC 310 is due to a runa- 
way effect. As particles accelerate to ultrahigh energies, electromagnetic 
cascades develop multiplying the number of charge carriers until their 
current shortcuts the gap. The excess particles are then swept away with 
the jet flow, until the gap reappears. 

Radio galaxies and blazars with very low accretion rates allow us to 
obtain a glimpse of the jet formation process near supermassive black 
holes. The sub-horizon variability in combination with the results from 
direct imaging campaigns invite to explore analogies with pulsars where 
particle acceleration takes place in two stages. In the first stage, particle 
acceleration occurs in the gaps of a charge-separated magnetosphere 
anchored in the ergosphere of a rotating black hole, and in a second stage 
at shock waves in the force-free wind beyond the outer light cylinder. 
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Fig. 1. The relation between Doppler factor and orientation 
angle evaluated for various values of the bulk Lorentz factor 
compared with the observational constraints of these 
parameters. The blue box shows the constraints on the 
Doppler factor arising from radio observations (32). For 
illustrative purposes, the red box shows the constraint from 
the gamma-ray optical depth to pair creation assuming 




Fig. 2. Significance map (color scale) of the 
Perseus cluster in gamma rays observed in 
the night of November 12/1 3th, 2012, with 
the MAGIC telescopes. The inset shows the 
radio jet image of IC 310 at 5.0 GHz obtained 
with the European VLB! Network (EVN) on 
October 29, 2012. Contour lines (and 
associated to them color scale) increase 
logarithmically by factors of 2 starting at 
three times the noise level (see supplement 
for image parameters). The ratio of the 
angular resolution between MAGIC and the 
EVN is 1:580 000. 
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Fig. 3. Average spectral energy distributions during the flare (red) along with 
previous measurements of IC 310 as observed by MAGIC. We show the results 
from the high (blue, open squares) and low (black, open markers) states reported 
in (28) and the average results (gray triangles) from (27) for the whole period. The 
dashed lines show power-law fits to the measured spectra, and the solid line with 
filled circles depicts the spectrum corrected for absorption in the extragalactic 
background light. As a reference, the spectral power-law fit of the Crab Nebula 
observations from (25) is shown (gray, solid line). Vertical error bars show 1 
standard deviation statistical uncertainity. Note that due to unfolding procedure 
spectral points are correlated. Horizontal error bars show the energy binning. 
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Fig. 4. Light curve^of IC 310 observed with the MAGIC telescopes in the night of 
November 12/13 , 201 2, above 300 GeV. As a flux reference, the two gray lines 
indicate levels of 1 and 5 times the flux level of the Crab Nebula, respectively. The 
precursor flare (MJD 56243.972-56243.994) has been fitted with a Gaussian 
distribution. Vertical error bars show 1 standard deviation statistical uncertainity. 
Horizontal error bars show the bin widths. 
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Fig. 5. Scenario for the magnetospheric origin of the 
gamma-rays: A maximally rotating black hole with event 

horizon (black sphere) accretes plasma from the center of 

the galaxy IC 310. In the apple-shaped ergosphere (blue) 

extending to 2r^ in the equatorial plane, Poynting flux is 

generated by the frame-dragging effect. The rotation of the 
black hole induces a charge-separated magnetosphere (red) 
with polar vacuum gap regions (yellow). In the gaps, the 
electric field of the magnetosphere has a component parallel 
to the magnetic field accelerating particles to ultra-relativistic 
energies. Inverse-Compton scattering and copious pair 
production due to interactions with low-energy thermal 
photons from the plasma accreted by the black hole leads to 
the observed gamma rays. 
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